Errors in particle tracking velocimetry with high-speed cameras 
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Velocity errors in particle tracking velocimetry (PTV) are studied. When using high-speed video 
cameras, the velocity error may increase at a high camera frame rate. This increase in velocity error 
is due to particle-position uncertainty, which is one of two sources of velocity errors studied here. 
The other source of error is particle acceleration, which has the opposite trend of diminishing at 
higher frame rates. Both kinds of errors can propagate into quantities calculated from velocity, such 
as the kinetic temperature of particles or correlation functions. As demonstrated in a dusty plasma 
experiment, the kinetic temperature of particles has no unique value when measured using PTV, but 
depends on the sampling time interval or frame rate. It is also shown that an artifact appears in an 
autocorrelation function computed from particle positions and velocities, and it becomes more severe 
when a small sampling-time interval is used. Schemes to reduce these errors are demonstrated. 

PACS numbers: 52.27.Lw, 52.70.Kz, 47.80.Cb, 07.05. Pj, 47.57.Gc 
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I. INTRODUCTION 

Particle tracking velocimetry (PTV) is a method to 
measure particle velocities fl| with video camera record- 
ing. In experiments, the particles are small solid objects 
that can scatter enough light to be imaged separately. 
The velocity is calculated based on measured positions 
of particles. Among the algorithms used to calculate ve- 
locities the most common is simply to divide the 
difference in a particle's position in two consecutive video 
frames by the time interval between the frames P, 0, @] . 
Although PTV may allow tracking individual particles 
for many frames, for this algorithm it is only necessary 
to track for two frames. 

Particle tracking velocimetry has been widely used for 
many years in topics in various fields, such as cell motion 
in biology Q, flow in granular materia ls B , and kinetic 
temperature in dusty plasmas 0, [9l-[l3l|7 Small parti- 
cles of solid matter can be added as tracers in a gas or 
liquid to study convection and turbulence [l] in fluid 
mechanics. Many dynamical quantities can be calculated 
using velocities measurements with PTV, for example ve- 
locity profiles in a shear flow mean-square vel ocity 
fluctuations [1], velocity distribution functions (9l.[l2l.ll5|. 
and velocity autocorrelation functions [3] ■ Some of these 
uses for PTV, such as velocity profiles, can also be accom- 
plished with particle image velocimetry (PIV) [13, [l8| . 
Compared to PIV, PTV provides a measurement of ve- 
locity at the location of a particle, without requiring an 
averaging over a grid. 

Recently, many experimenters using PTV have taken 
advantage of the abundance of high-speed cameras now 
offered for commercial sale. They may be unaware, how- 
ever, that velocities determined using the PTV method 
can have errors that become more severe as the cam- 
era frame rate is increased. For example, we show that 
when determining kinetic temperature by calculating the 
mean-square velocity fluctuation for random motion, the 



result will have an exaggerated value that worsens at 
higher frame rates, due to one of the two kinds of er- 
rors studied here. 



II. VELOCITY ERRORS IN PTV 

We identify two kinds of velocity errors in PTV accord- 
ing to their source. One kind arises from acceleration of 
the particle during the time interval between measure- 
ments of its position. The other arises from errors in the 
position measurements themselves. While the former is 
made less severe by using a faster frame rate, the latter 
is actually made worse. Therefore, a faster frame rate is 
not always best for PTV. 



A. Velocity error arising from acceleration 

An error arises from acceleration, whether due to a 
change in a particle's speed or its direction. This error 
occurs for all kinds of acceleration, for example a particle 
accelerating along a straight line, moving in a circle at a 
steady speed, or colliding with another particle. 

It is impossible to eliminate this error when the motion 
is unknown during the sampling interval between frames. 
The only information available in PTV are the particle 
positions determined at the times that a video image was 
recorded, i.e., the video frames. 

The size of this error also depends on the algorithm 
for calculating velocity from data for the particle posi- 
tions. The simplest and most common algorithm is to 
assume that the particle moves with zero acceleration be- 
tween two consecutive frames. In this two-frame tracking 
method, the particle is assumed to move in a straight line 
at a steady speed between the positions in two frames. 
Thus, the velocity is calculated simply as the difference 
in positions divided by the time interval between frames. 
Using two position measurements Xmcas, denoted by j 
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and j + 1 separated by a time At, the velocity at a time 
halfway between them is calculated as 

The camera's frame rate is 1/Ai if no frames are skipped. 
We use upper-case symbols like V and X to indicate 
quantities that are computed from the positions that are 
measured in a single frame. If the particle actually un- 
dergoes acceleration, for example if it has a curved trajec- 
tory, the simple algorithm of Eq. (IT|) will obviously lead to 
errors in the velocity V. Alternatively, one could use an 
algorithm using more than two position measurements, 
which can sometimes better account for acceleration. A 
spline fit could reduce the error arising from accel- 
eration, but it is more computationally expensive than 
two-frame tracking. 

To illustrate how errors in V arise from acceleration 
and how these errors diminish with a higher frame rate, 
we consider the motion during a Coulomb collision of 
a pair of identical electrically-charged particles. The 
largest acceleration in this case will occur when the par- 
ticles are closest and the direction of motion is changing 
most rapidly. This highly curved portion of a trajectory 
can be approximated as a circular arc. This motivates 
us to choose uniform circular motion as a simple instruc- 
tive example of velocity errors arising from acceleration. 
Fig. 1. Suppose the particle's position is measured accu- 
rately at each frame, as indicated by the dots in Fig. 1(a). 
The simple two-frame tracking method, which assumes 
zero acceleration between measurements, will describe 
the motion as a polygon instead of the ideal circle. At 
a fast frame rate, the polygon has more sides and more 
closely approximates the ideal circle. 



B. Velocity error arising from position uncertainty 

Another error in velocity arises from uncertainties in 
the particle positions from which the velocity is com- 
puted. Particle positions are uncertain for at least two 
reasons: random noise in the camera sensor, and the fi- 
nite size of pixels in the sensor. The latter leads to the 
phenomenon of pixel locking, where particles are wrongly 
identified as being located at favored positions such as 
the corner or middle of pixels. While it is possible to 
design an experiment to reduce these particle-position 
uncertainties [J, [ist , they can never be eliminated. 

Suppose that all measured particle positions contain 
an uncertainty Sx. We can use propagation of errors 
to estimate the uncertainty of the calculated velocity V, 
arising from the uncertainties in the particle positions, 
for any given algorithm. Choosing the algorithm used 
in the simple two- frame tracking method, Eq. (IT|), the 



uncertainty in the calculated velocity V is 

Importantly, this source of error becomes larger, not 
smaller, as the time interval between measurements At 
is decreased. This is seen in Eq. where the denom- 
inator diminishes with an increasing frame rate but the 
numerator does not vary with frame rate at all. Thus, 
this error in velocity is cx 1/At. Now that high frame- 
rate scientific cameras have become more easily available, 
this error in PTV probably occurs more commonly. 

To illustrate the combined effect of both types of er- 
rors, those due to acceleration and due to uncertainties 
in particle position, consider the sketch in Fig. 1(b) for 
a particle undergoing uniform circular motion. The un- 
certainty in particle positions is indicated in Fig. 1(b) 
by shading around the true positions. Uncertainties in 
particle position cause the measurements to fall on the 
vertices of a polygon that is deformed, as compared to 
Fig. 1(a) without particle-position uncertainties. Here we 
are interested primarily in errors in the velocity V. In 
Fig. 1(b), the length of the edges of the distorted poly- 
gon is an indication of velocity. Comparing this thick 
irregular polygon to the original ideal circle, it is obvious 
that the velocity vector will have errors in its direction. 
The magnitude of the velocity will also have errors, as 
we discuss next. 

C. Illustration of PTV velocity errors for a particle 
in uniform circular motion 

To demonstrate the combined effect of both sources of 
error in velocity, we present a simple simulation. A single 
particle is assumed to perform uniform circular motion 
with a radius R and period T. A time series of a parti- 
cle's X and y coordinates is recorded at intervals At. The 
simulation duration is 1000 circular periods. To simulate 
a measurement error, we add a random error chosen from 
a Gaussian distribution to the true position. The prob- 
ability of an error in this Gaussian distribution is 
oc ex-p[—xl^^/{2Sl)], and the same for the y direction. We 
then calculate a time series V(t), for the particle's veloc- 
ity using the simple two-frame tracking method, Eq. ([l}. 
We use V(t) to calculate a time series for the kinetic en- 
ergy, KE{t) = TO|V(t)p/2, where the mass m will cancel 
in our final results. We average over the entire time se- 
ries yielding (KE), which we compare to the true kinetic 
energy KEtruc in Fig. 2. We indicate the discrepancy be- 
tween (KE) and KEtruc as a single data point in Fig. 2. 
We vary the size of the position error and the sampling 
interval At, yielding all the data points shown in Fig. 2. 

Figure 2 shows the total error in (KE) due to both 
sources combined. For 6x = 0, Fig. 2 also shows the er- 
ror due to acceleration only, as indicated by circles. This 
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error due to acceleration is always negative, and it is most 
severe for large At, i.e., for a slow frame rate. The total 
error in (KE), however, can be either positive or negative, 
as shown by the other data points in Fig. 2. The total er- 
ror is positive at small At due mainly to the contribution 
of particle-position uncertainty, and it becomes negative 
at large At due mainly to acceleration. In between, the 
total error in (KE) has its smallest magnitude, which for 
this simulation occurs at at a sampling interval of about 
3% of the circular period, i.e., At/T w 0.03. This obser- 
vation suggests that it may be possible to choose a best 
frame rate to minimize errors, as we discuss next. 

In an actual experiment, one cannot independently 
measure the two sources of error in velocity, or in quan- 
tities such as (KE) that are calculated from the velocity. 
In most physical systems, the true motion of one particle 
will differ from that of another particle, and their accel- 
erations will vary with time, unlike the idealized circu- 
lar motion simulation in Fig. 2. Nevertheless, examining 
Fig. 2 suggests a possible scheme for choosing a At that 
has some promise to reduce the total error. 



D. Scheme for reducing total error 

The scheme we suggest for reducing the total error 
for a quantity (such as (KE)) computed from velocities 
is as follows. The experimenter can record motion at a 
high frame rate, and then analyze the data not only at 
the frame rate, but also at larger At by skipping one 
frame to double At, two frames to triple At, and so on. 
Plotting the average of KE vs. sampling interval At will 
yield a graph similar to Fig. 2. If the graph exhibits a 
nearly flat spot at a particular value of At, that value is 
a candidate to consider as the best value for computing 
KE. In Sec. Ill B, we will examine actual experimental 
data as a demonstration. 



III. PTV VELOCITY ERRORS IN 
EXPERIMENTS 

A. Experiment 

Various experiments using PTV can have much in com- 
mon, even when the physical systems being studied are 
completely different. Particle positions are measured in 
images that correspond to individual video frames, and 
to do this they must be spaced sufficiently that they can 
be distinguished. In a granular media experiment, the 
particles might fill about half the volume 1] , while in a 
dusty plasma experiment they are more widely spaced, 
filling typically one millionth of the total volume. There 
are several choices for image analysis methods to mea- 
sure the particle position, which have various advantages 
for different physical systems; for our dusty plasma we 



use the moment method of What might differ the 
most, from one physical system to another, is the na- 
ture of the true particle motion. In granular media, par- 
ticles interact as hard spheres, with large accelerations 
during the brief collision events. Dust tracers in fluid 
mechanics experiments represent the opposite extreme: 
they do not collide at all, but are merely swept along 
with the flow of the fluid. Between these two extremes 
is dusty plasma, with particles that collide softly with 
a screened Coulomb repulsion, so that they undergo a 
gradually changing acceleration. For recording videos of 
physical systems where the particles do collide, it is use- 
ful to characterize a typical time scale. A typical time 
scale is w 10~^ s between hard-sphere collisions in granu- 
lar flow Q . The duration of a Coulomb collision in dusty 
plasma is « 10"^ s jT9||. 

We carry out a demonstration experiment [20j us- 
ing PTV. Polymer microspheres (which we will refer to 
as dust particles) of 8.1-/xm diameter are immersed in 
a partially- ionized argon gas under vacuum conditions. 
This four-component mixture (micron-size dust parti- 
cles, rarefied gas atoms, electrons and positive ions) is 
called a dusty plasma. The dust particles are electri- 
cally charged, and they are levitated against the down- 
ward force of gravity by a vertical electric field in the 
plasma. Fig. 3(a,b). Due to their mutual Coulomb re- 
pulsion, the dust particles tend to arrange themselves so 
that they are separated by a distance b w 0.7 mm, which 
is much l arg er than their diameter. Using a laser heating 



method [15|, |20| , the kinetic energy of dust particles is 



augmented, so that their velocities fluctuate with a typi- 
cal magnitude of « 1 mm/s, as we will find in Sec. Ill B. 
The dust particles always remain within a single horizon- 
tal layer, which is well suited for imaging. 

This horizontal layer of dust particles is imaged by a 
high-speed camera viewing from above. We use a 12-bit 
Phantom v5.2 camera, with 1152 x 720 pixel resolution, 
and a Nikon 105 mm focal-length macro lens. The dust 
particles are illuminated by a horizontal sheet of 576 nm 
laser light, and the camera lens is fitted with a filter to 
block unwanted light at other wavelengths. The camera's 
36.2 X 22.6 mm2 field-of-view (FOV) includes N « 2100 
dust particles. A video is recorded as a series of bit-map 
images. The image of one dust particle fills about ten 
pixels. Fig. 3(c). Due to the large interparticle spacing, 
only about 2% of the pixels in an image are brighter than 
the background. The camera recorded a 20 s movie at 
250 frames/s, i.e., a time between frames of 4 ms. 

We can compare our 4 ms time interval between frames 
to two physical time scales in the experiment, which are 
both of interest in determining the velocity error. One 
physical time scale is for acceleration for a Coulomb col- 
lision amongst dust particles. This is typically « 10"^ s, 
based on time required for the dust-particle velocity to 
change significantly in the binary-collision experiment of 
Konopka et al. |l9|. The other physical time scale (for 
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errors arising from particle-position uncertainty) is the 
time required for a dust particle to move one pixel. A 
thermal velocity of w 1 mm/s is a typical velocity in 
our experiment, and at this speed a dust particle would 
require w 30 ms to cross one pixel. 

Images for each frame in the movie are analyzed us- 
ing the method of [ij to measure the dust particle posi- 
tions. We compute velocities using the simple two-frame 
tracking method, Eq. ([T]). We then calculate the kinetic 
temperature and a correlation function, as we describe 
next. 



B. Kinetic temperature 

We will demonstrate here, using experimental data, 
that PTV does not yield a unique value for the kinetic 
temperature. Instead, it yields a value that depends on 
the sampling interval At between images that are ana- 
lyzed. All the data presented below come from the same 
experiment, with the camera always operated at the same 
frame rate. To demonstrate the dependence on the sam- 
pling interval, we will repeat our analysis by skipping 
frames. 

The quantity that we test here is the kinetic temper- 
ature, averaged over a time series. The time series for 
kinetic temperature is calculated from mean-square ve- 
locity fluctuations 
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where N is the number of dust particles analyzed. (Al- 
ternatively, one could calculate the temperature as a 
fit parameter assuming a Maxwellian velocity distribu- 
tion [l^.) We then average over the entire time series, 
yielding the calculated result for the kinetic temperature 
(T). In Eq. ([3]), V(t) indicates a velocity averaged over 
all dust particles in one frame, i.e., the center-of-mass 
velocity. (The bar indicates an average over all dust par- 
ticles in the FOV in one image, while the brackets (• • •) 
indicate averages over time.) 

The results for the calculated kinetic temperature, as 
the sampling time interval At is varied, are shown in 
Fig. 4. We see that (T) decreases gradually as the sam- 
pling time interval increases. The error in kinetic tem- 
perature is the difference of (T), which varies with At, 
and an unknown true value, which does not vary. Thus, 
it is clear that PTV does not yield a unique value for the 
kinetic temperature, but instead a value with an error 
that depends on the experimenter's choice of At. More- 
over, the error is not necessarily reduced by choosing a 
faster frame rate and using every frame without skipping 
frames; if one did this, the result for (T) would become 
steadily larger as the frame rate is increased. 

Having observed experimentally that there must be an 
error in (T) that depends on the experimenter's choice 



of At, we ask how the experimenter should choose At. 
Motivated by Sec. II D, we will use our flat-spot scheme 
with our graph of (T) vs. At, Fig. 4. 

Examining Fig. 4, we see a general trend of (T) di- 
minishing with the sampling time interval At. For the 
smallest At shown, the curve has a steep slope, while 
for larger At there is what appears to be a nearly flat 
spot. Comparing to Fig. 2 for ideal circular motion, we 
interpret the steep slope at small At as an indication 
that the error is dominated by particle-position uncer- 
tainties, which is the source of error that becomes worse 
with higher frame rates. In the nearly flat spot of Fig. 4 
for our experiment, we would choose At ~ 0.03 s, corre- 
sponding to (r) « 2.5 X 10^ K and thermal velocity is 
1.3 mm/s. 

A limitation of this scheme for choosing At is that 
the dust particles in the experiment do not undergo the 
same acceleration all the time, unlike the simple circular- 
motion that is simulated in Sec. II C. Therefore, there is 
no strong reason to expect that the nearly flat spot in 
the graph of (T) vs. At will be the same as in Fig. 2 for 
the circular-motion simulation. Nevertheless, this seems 
to be the best guidance available to us, given what we 
know about the errors in V and therefore (T) that arise 
from particle-position uncertainty and acceleration. 

In addition to the flat-spot scheme described above, 
we can suggest an alternate scheme which leads to the 
same choice of At in the case of our experiment. As de- 
scribed in Sec. Ill A, there are two physical time scales 
that can be compared to At. As an upper limit for At, to 
avoid excessive velocity errors arising from acceleration, 
it is desirable to choose At significantly smaller than the 
time scale for velocity to change significantly in the phys- 
ical system. As a lower limit for At, to avoid large errors 
arising from particle-position uncertainty. At should be 
comparable to, or larger than, the time required for a 
particle at a typical velocity (in our experiment the ther- 
mal velocity) to move one pixel. Our reasoning behind 
the lower limit is that the particle-position uncertainty 
is generally a certain fraction of one pixel [l|, [l^l , so that 
if the displacement of a particle in At is much smaller 
than one pixel a large velocity error arising from particle- 
position uncertainty should be expected. For this ex- 
periment, the upper limit due to acceleration is about 
100 ms, while the lower limit due to particle-position 
uncertainty is about 30 ms. Choosing At significantly 
smaller than the upper limit, and comparable to the lower 
limit, leads us to choose the same value as the flat-spot 
scheme. At w 0.03 s in our experiment. 

We note that another approach to reduce the effect of 
particle-position uncertainty was reported by Knapek et 
al. 21|. In a dusty plasma experiment similar to ours. 



they operated a camera at 500 frames/s. Before carrying 
out computations to track the particles and compute the 
kinetic temperature, they superimposed three consecu- 
tive bit-map images, averaging them pixel by pixel. This 
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averaging of images has the effect of reducing particle- 
position uncertainty, with the tradeoff of a reduction of 
time resolution. In our scheme, instead of averaging in- 
formation while it is still in the form of bit-map images, 
we focus our efforts on the particle-tracking algorithm 
and choosing the time interval. 



C. Transverse current autocorrelation function 

Aside from the kinetic energy and related quantities 
like the kinetic temperature, there are other quantities 
one might wish to compute from velocity data produced 
by PTV. For a correlation function computed from exper- 
imental velocity and position data, we will demonstrate 
here that artifacts can arise at high frame rates that are 
different from the errors that appear in calculations of the 
kinetic temperature. We will also demonstrate that re- 
ducing these artifacts can require a scheme different from 
the one described in Sec. II D for kinetic temperature. 

The correlation function we will choose is the trans- 
verse current autocorrelation function Cxit), which is 
used in the study of shear motion and viscoelastic- 
ity [iO, 22, [2^. This correlation function is computed 



from records of particle positions Xi{t) and velocities 
Vy^i{t) in orthogonal directions. 



Crit) = {j;{k,o)jy{k,t))/{j;{k,o)jy{k,o)), (4) 



where 



EN 
Vy i{t) eyip[ikxi{t)] 
i—l 



(5) 



The quantity jy{k,t) is called a "transverse current," 
even though it has nothing to do with electric current. 
Both the current and the correlation function Crit) de- 
pend on a parameter k that is adjustable. In the study 
of viscoelasticity [13], values of k comparable to 3/6 are 
typical, where 6 is a mean interparticle distance. (In our 
experiment, b ~ 0.7 mm.) Like many autocorrelation 
functions, Crit) has an initial decay at small time t; this 
initial decay is of great interest in the study of viscoelas- 
ticity [20]. 

Using particle velocity and position data from our ex- 
periment, we find an undesired artifact in the initial 
decay of the calculated Crit), which can be seen in 
Fig. 5(a). This artifact is most prominent when the 
camera frame rate is high and when we use the simple 
two-frame tracking method to prepare the velocity data 
used in calculating Crit). We attribute this artifact to 
particle-position uncertainty. 

Although this artifact arises due to a high camera 
frame rate, it is unattractive to attempt to eliminate it 
by reducing the frame rate. Many data points during the 
initial decay are required for the study of viscoelasticity, 
and this requires a high frame rate. Thus, we need a 



different approach to reduce the artifact, like the two- 
step approach we describe next, for reducing the effects 
of particle-position uncertainty. 

In the first step, we use a three-frame tracking method, 
which reduces the artifact as seen in Fig. 5(b). In this 
three-frame tracking, we do not entirely skip any frames, 
because that would reduce the temporal resolution of 
Crit). Instead, we interlace pairs of frames, and cal- 
culate the positions and velocities as: 



= (xj_i + + Xj+i)/3, 



V, = (x,+i -x,_i)/2Ai. 



(6) 



(7) 



Here, the positions of a particle in any 

three consecutive frames. Equations ^ and ([7]) can be 
derived from the linear regression for three data points. 
Using three frames rather than two requires a sufficiently 
fast frame rate, which is satisfied in our experiment with 
a time interval between frames of 4 ms, which is much 
shorter than the time scale for changes in the dust parti- 
cle's velocity w 10~^ s. (While this three-frame tracking 
method has the advantage of reducing the unwanted ar- 
tifact in the initial decay of Crit), it may also have an 
effect on the long-time oscillations of Crit), which we 
have not studied.) 

To further reduce the artifact, in the second step, we 
smooth jy{k,t) before calculating Crit). This averaging 
requires a sufficiently high frame rate. We smooth the 
transverse current jy{k,t) with a moving average over 
five consecutive frames before calculating CT{t). The 
result for Crit), Fig. 5(c), is a great reduction of the 
artifact without an unnecessary loss of time resolution in 
the initial decay. 

To verify that these two steps reliably remove artifacts 
without changing the true shape of the initial decay, we 
carried out numerical simulations, where the true shape 
is known. These simulations used a Langevin molecu- 
lar dynamic method (20l | to track particles by integrating 
their equations of motion. To mimic the experiment's 
particle-position uncertainty, we added random errors to 
the positions in the simulation as in Sec. II C. The sim- 
ulation results, not shown here, confirm that the artifact 
is diminished by the two steps described above, and they 
also confirm that this improvement results in only a neg- 
ligible change in the initial decay of the calculated Crit), 
as compared to the true Cxit) without the random er- 
rors. 



IV. SUMMARY 

In summary, we have studied two sources of errors in 
PTV using high-speed cameras. What may surprise some 
users is that a higher frame rate is not necessarily better, 
and in fact it can greatly worsen velocity errors arising 
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from particle-position uncertainty. We described our so- 
lutions to reduce these errors: two schemes that lead to 
the same choice for Ai when calculating kinetic tempera- 
ture, and a two-step approach for computing an autocor- 
relation function. We demonstrated these solutions using 
the data from a dusty plasma experiment. This work was 
supported by NSF and NASA. 
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FIG. 1: (Color online) Simulations of particle trajectories for 
circular motion (a) with no particle-position uncertainty, and 
(b) with particle-position uncertainty. The true particle po- 
sitions at time intervals At are marked by circular dots. The 
triangles in (b) indicates particle positions including an error 
from a Gaussian distribution with an uncertainty S^- 



FIG. 2: (Color online) Velocity errors for a simulated par- 
ticle executing uniform circular motion with period T, as 
in Fig. 1(b). The velocity is calculated using Eq. ([T]), and 
then the kinetic energy is calculated and averaged over 1000 
cycles of the circular motion. The true kinetic energy is 
m{2-KR/Tf /2. The severity of the two sources of error in 
the kinetic energy depends on the sampling interval At. The 
dominant error arises at small At from particle-position un- 
certainty, and at large At from acceleration. There is a nearly 
flat spot at At/T ~ 0.3, where the total error from two sources 
is smallest. 



FIG. 3: (Color online) Configuration for the experiment, (a) 
Dust particles were levitated in a horizontal layer above a 
lower electrode. A vacuum chamber, not shown here, was 
filled by a weakly-ionized plasma, consisting of electrons, pos- 
itive argon ions, and a neutral argon gas. A pair of 532-nm 
laser beams accelerated the dust particles to raise the kinetic 
temperature. A top-view camera imaged the dust particles, 
which were illuminated by a horizontal sheet of 576-nm laser 
light, (b) Dust particles are electrically levitated in a hori- 
zontal layer, (c) Experimental bitmap image from one frame 
of a video. Only ~ 1/60 of the experimental image is shown. 
The inset is a magnified image of one dust particle. 



FIG. 4: Experimental data for the kinetic temperature, calcu- 
lated using Eq. Q. To vary the sampling time interval At, we 
skipped frames in the data analysis. We note the appearance 
of a nearly flat spot, similar to the one seen in Fig. 2 for the 
circular-motion simulation. In Sec. Ill B we suggest that for 
calculating kinetic temperature, it is desirable either to choose 
At in the flat spot of this graph, or to choose At comparable 
to the time required for a particle at the thermal velocity to 
move a distance of one pixel. Both of these schemes lead to 
the same choice in this experiment, At « 0.03 s. 



FIG. 5: Experimental data for Crit), an autocorrelation func- 
tion of the time series jy{k,t), as defined in Eq. An un- 
wanted artifact appears in the initial decay in (a) with two- 
frame tracking, i.e., using Eq. ([T}) to calculate velocities. This 
artifact is reduced in (b) by using three-frame tracking, i.e., 
using Eqs. ((Sjl and ([7]). It is further reduced in (c) by smooth- 
ing the the time series jy {k, t) before calculating its autocor- 
relation. Data shown here were computed for kb = 3.5, where 
h is the mean interparticle distance. 



